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O : ABSTRACT 

Different wavelength, regimes and methods for estimating the space density of the star- 
formation rate (SFR) result in discrepant values. While it is recognised that ultra-violet (UV) 
^ , and Ha emission line data must be corrected for the effects of extinction, the magnitude of the 

required correction is uncertain. Even when these corrections are made there remains a signif- 
icant discrepancy between SFRs derived from UV and Ha measurements compared with those 
derived from far-infrared (FIR) and radio luminosities. Since the FIR/radio derived SFRs are 
^3 . not affected by extinction, and simple corrections to reconcile the UV and Ha measurement 

with these do not fully account for the discrepancies, a more sophisticated correction may be 
required. Recent results suggest that at least part of the solution may be a form of extinction 
' which increases with increasing SFR (or luminosity, given the common assumption that SFR is 

proportional to luminosity). We present an analysis of the effects of a dust reddening dependent 
on star formation rate applied to estimators of SFR. We show (1) that the discrepancies between 
Ha and FIR/radio SFR estimates may be explained by such an effect, and we present an iterative 
method for applying the correction; and (2) UV-based estimates of SFR are harder to reconcile 
with FIR/radio estimates using this method, although the extent of the remaining discrepancy is 
less than for a non-SFR-dependent correction. Particularly at high redshift, our understanding 
of extinction at UV wavelengths may require a still more complex explanation. 

Subject headings: dust, extinction — galaxies: evolution — galaxies: starburst 

1. Introduction history of the universe. It is generally agreed that 

p* , the comoving space density of the star forma- 
In recent years a variety of observations and ob- ^ion rate (SFR) in galaxies, rises by an order of 

servational methods have been used to contribute magnitude between ^ = and z = 1 (Lilly et 

to our understanding of the global star formation ^1. 1996; Connolly et al. 1997; Hogg et al. 1998; 

; Flores et al. 1999; Haarsma et al. 2000). Despite 
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z > 2 however, (Madau et al. 1996; Connolly et 
al. 1997; Hughes et al. 1998; Yan et al. 1999; Stei- 
del et al. 1999; Hopkins et al. 2000) it is still un- 
clear whether the evolution of p* reaches a peak 
around 2 « 1.5 and decreases significantly there- 
after (e.g., Madau et al. 1996), or whether it stays 
flat to much higher redshifts (e.g., Steidel et al. 
1999; Haarsma et al. 2000). 

Some of this uncertainty may be related to the 
difficulty of making corrections for extinction at 
ultra-violet (UV), optical and near-infrared wave- 
lengths. In the absence of extinction corrections 
there is a clear discrepancy between values of p* es- 
timated from observations at radio or far-infrared 
(FIR) wavelengths and those estimated from Ha 
and UV measurements. Even when simple extinc- 
tion corrections to the UV and Ha measurements 
are made, however, systematic differences persist 
between estimates of SFR or p* made at differ- 
ing wavelengths. This can be seen, for example, 
in the distribution of SFRs for (1) individual ob- 
jects in a study of local galaxies by Cram et al. 
(1998, their Figure 1), and in a mult i wavelength 
study of Selected Area 57 (SA57) by Sullivan et 
al. (2001, their Figure 3); and (2) globally, in a 
diagram showing the redshift dependence of p* by 
Haarsma et al. (2000, their Figure 10). 

Several efforts have been made to address this 
concern. A detailed investigation of some of these 
issues is presented in Adelberger & Steidel (2000), 
with particular emphasis on obscuration at UV 
wavelengths. Calzetti & Heckman (1999) have 
modelled the evolution with redshift of galaxy 
opacity for applying corrections to UV-derived es- 
timates of p* . Meurer et al. (1999) have refined 
a technique for correcting UV luminosities for 
attenuation using rest-frame luminosities alone. 
This method is based on an observed correla- 
tion between the FIR/UV flux ratio and j3 the 
UV spectral slope. Gordon et al. (2000) presents 
a more general flux ratio method for establish- 
ing wavelength-specific attenuations for individual 
galaxies, based on assumed SEDs. 

We extend these investigations by first deriv- 
ing an SFR-dependent (or luminosity-dependent) 
attenuation from existing observations. We then 
examine its efficacy by applying it to Ha and UV 
estimates of SFR for a large sample of local galax- 
ies. The examination is subsequently expanded to 
investigate the effects of the application to global 



SFR densities, derived from Ha or UV measure- 
ments, spanning a broad redshift range. 

In Section 2 we comment on the use of FIR 
and 1.4 GHz luminosities as estimators of SFR. 
Section 3 describes the formulation of an SFR- 
dependent reddening parametrisation. This is 
then applied to correct estimates of SFR for a sam- 
ple of local galaxies in Section 4, and subsequently 
to estimates of p* spanning a broad redshift range 
in Section 5. Scc;tion 6 presents a brief investiga- 
tion of the empirical relation between SFRuv and 
intrinsic SFR, and its effects when used blindly to 
correct estimates of at all redshifts. 

Throughout this paper we have used SFR cal- 
ibrations defined consistently to account for the 
mass range 0.1 < M < 100 assuming a 
Salpeter IMF. We adopt the Ha and UV SFR 
calibrations given by Kennicutt (1998), shown in 
Table 1. Our assumed FIR and 1.4 GHz SFR cal- 
ibrations are also given in this Table, and are ex- 
plained in more detail below. We assume = 
75kms~^Mpc~^ and go = 0.5, and have converted 
previously published estimates of SFR to our as- 
sumed calibrations and this cosmology, where nec- 
essary, to ensure consistency. 

2. FIR and Radio estimates of SFR 

The intrinsic SFR of a galaxy cannot be esti- 
mated directly from observations at UV or optical 
wavelengths because of the uncertainties caused 
by the presence of unknown amounts of obscuring 
dust. Indicators of SFR at longer wavelengths, 
at which dust is transparent, are sought as an 
alternative. In the past few decades FIR lumi- 
nosity has been used extensively as an estima- 
tor of current star formation rate in galaxies, and 
throughout this paper we have adopted the cali- 
bration given by Kennicutt (1998), which is shown 
in Table 1. While FIR luminosity is expected to 
be an excellent SFR tracer for strong, compact, 
dusty starbursts, the situation is less clear when 
the disks of normal galaxies are being examined. 
The uncertainties in using FIR luminosity as an 
SFR estimator for normal galaxies are described 
by Kennicutt (1998), and while FIR luminosity in 
late-type star-forming galaxies appears to corre- 
late well with other estimators of SFR, this is not 
generally the case for early-type spirals and ellip- 
ticals. 
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Table 1 

Assumed luminosity-to-SFR calibrations. 



Indicator Calibration 



UV SFRuv = iuv/7.14xl020WHz-i 

Ha SFRHa = iHa/1-27 X 1034 W 

FIR SFRfir = Lfir/5.81 x IO^Lq 

1.4 GHz SFRi.4GHz = Li.4GHz/8.4 x lO^o W Hz"^ 



At still longer wavelengths, populations of star- 
forming galaxies have been detected in increas- 
ingly sensitive radio surveys (e.g. Windhorst et al. 
1984, 1985; Oort 1987; Hopkins ct al. 1999). As 
well as supporting the well-established radio-FIR 
correlation for disk or star-forming galaxies (which 
extends over more than 4 orders of magnitude and 
is commonly accepted to be a product of star for- 
mation processes, see review by Condon 1992), es- 
timates of SFR from measurements of 1.4 GHz ra- 
dio luminosity are now becoming more common 
(e.g. Condon 1992; Cram et al. 1998; Haarsma 
et al. 2000). We have adopted the 1.4 GHz SFR 
calibration as given by Haarsma et al. (2000), 
SFRi.4GHz = Q(Li.4/4.6 X lO^iWHz"^), with 
Q = 5.5 the factor accounting for the inclusion 
of stars less massive than 5 Mq (sec Tabic 1). 

As a result, the use of radio or FIR luminosi- 
ties as estimators of SFR appear to be justified for 
samples of galaxic^s which lie on the radio-FIR cor- 
relation. Admittedly, there are still concerns over 
the details of the physical processes which pro- 
duce this correlation (Condon 1992). While the 
currently favoured explanation are processes re- 
lated to star formation, some radio-quiet quasars 
also appear to follow the correlation (Cram et al. 
1992), as do Seyfert galaxies lacking compact ra- 
dio cores (Roy et al. 1998, although with greater 
scatter in this case). The scatter in the radio-FIR 
correlation for disk galaxies, while smaller, is also 
significant and its origins remain unknown. De- 
spite these concerns, long-wavelength estimators 
of SFR provide a valuable tool which avoids the 
effects of obscuration, thus allowing useful insights 
into the properties of star formation in various 
galaxy populations. 



3. An SFR-dependent reddening 

Star forming galaxies in the local universe have 
been shown to exhibit a correlation between ob- 
scuration and FIR luminosity, as seen for example 
in Wang & Heckman (1996). This implies that 
the obscuration in a galaxy is related to its SFR. 
This same general effect can also be seen in other 
recent results, e.g., Adelberger & Steidel (2000), 
particularly their Figure 11, Calzetti et al. (1995), 
and Sullivan et al. (2001). 

Wang & Heckman (1996) demonstrate a trend 
between FIR luminosity (Lfir) and the inverse of 
the Balmer decrement for a small population of 
normal Hubble types as well as some Markarian 
starbursts. A similarly small sample of nuclear 
starburst and blue compact galaxies analysed by 
Calzetti et al. (1995) show a clear trend between 
colour excess, (or UV spectral index), and -Lfir- 
Both trends are in the sense of increasing obscu- 
ration with increasing ipiR- The scatter in these 
trends is notable, however, and may be related to 
intrinsic differences within galaxy populations, as 
well as between the different populations analysed. 

Now, since FIR luminosity for these objects can 
be treated as an estimator of a galaxy's intrinsic 
star formation rate, this suggests that these rela- 
tions can be recast in terms of obscuration (colour 
excess, for example) and SFRfir. Using the data 
of Wang & Heckman (1996) we perform a least 
squares fit to the Balmer decrement (not its in- 
verse) and log(LFiR) to derive 



0.797 log 



FIR 



3.952. 



(1) 



Fh0 \ Lq 

The Pearson correlation coefficient for this fit is 
0.6. Sullivan et al. (2001) independently find a 
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relation consistent with this result, from observed 
Balmer decrements for their UV (200 nm) selected- 
sample. E{B — V), the colour excess appropriate 
for nebular emission lines, is then calculated using 
the Balmer decrement (Calzetti et al. 1996)'^: 



E{B - = 



log(i?obs/i?int) 



(2) 



Here i?obs and iiint(=2.88 for Case B recombina- 
tion) are the observed and intrinsic values for the 
Balmer decrement. We can combine Equations 1 
and 2 to form a direct relationship between colour 
excess and FIR luminosity. Figure 1 shows the 
resulting relation, along with the data used to de- 
rive Equation 1, (modified from Figure 10 of Wang 
& Heckman 1996), and data taken from the in- 
dependent sample of starburst and blue compact 
galaxies by Calzetti et al. (1995). Since a negative 
colour excess is unphysical, we assume that lumi- 
nosities which give rise to E(B — V) < from this 
relation correspond to zero attenuation, and this 
is seen in the flattening of the line in Figure 1 . Be- 
low LpiR ~ 4 X 10"'^'' Lq the data from Calzetti et 
al. (1995) appear to be consistent with the trend 
seen in the Wang & Heckman (1996) data, but 
at higher luminosities even greater colour excesses 
are seen. This suggests that for high-SFR galaxies 
in the local universe, the following analysis may be 
somewhat conservative. 

The relation between colour excess and FIR lu- 
minosity can now be transformed to one between 
colour excess and SFR, using our adopted cali- 
bration for SFRpiR, and we then find the colour 
excess appropriate to a given intrinsic SFR: 

0.797 log(SFRFiR) 3.834 



£;(B-y)gas = 1.965 log 

z.6o 

(3) 

This, then, forms the crux of an SFR-dependent 
reddening relation. Given the intrinsic SFR for a 
galaxy (which may be derived from radio as well 



^We use the reddening curve k{\) = A{X} / E{B — V) for 
star-forming systems formulated by Calzetti et al. (2000): 

fe(A) = 2.659(-l. 857 + 1.040/A) -1-4.05 

(0.63 /im < A < 2.20 /im) 
= 2.659(-2.156 + 1.509/A - O.IQS/A^ 
-1-0.01 1/A^) -1-4.05 

(0.12 /im < A < 0.63 /im) 



as FIR luminosities, for example), and a reddening 
curve, the attenuation to be applied at all wave- 
lengths can be calculated. Using the standard for- 
mulation (e.g. Calzetti et al. 2000; Calzetti 1997) 
the intrinsic luminosity, Li{\), is related to the 
observed value, Lo{\), by: 



= io(A)100-4^(^- 



y)fc(A) 



(4) 



and since all the assumed SFR calibrations are di- 
rectly proportional to luminosity, this can be con- 
sidered to be a relation between intrinsic and ob- 
served SFR. 

Additionally, and most importantly, since 
Equation 3 is monotonic even a measurement of 
the attenuated SFR can be used to estimate the 
amount of reddening. This is done by substitut- 
ing Equation 3 into Equation 4 (cast in terms of 
SFR rather than luminosity, since the SFR cali- 
bration constant cancels) and solving numerically 
the resulting transcendental equation, which for 
the wavelength of Ha can be written: 



\og{SFRi) = 



log(5Fi?„(Ha)) + 2.614 X 

, r0.7971og(SFRi) 3.834 
log 



(5) 



Hence, by assuming such a form for an SFR- 
dependent obscuration, even observations of a sin- 
gle Ha emission line (or other star-formation- 
sensitive line, such as [Oil]) may be used to es- 
timate attenuation and intrinsic SFR. 

Now the colour excess appropriate for the stel- 
lar continuum, E{B — V)star, is related to that for 
nebular gas emission lines by 



E{B - y),t,, = O.UE{B - V) 



gas; 



(6) 



(e.g. Calzetti et al. 2000; Calzetti 1997). As a re- 
sult it should in principle be possible to apply the 
same procedure to derive appropriate corrections 
for UV-continuum estimates of SFR. Combining 
Equations 6, 3 and 4 gives the following relation 
for UV wavelengths: 



\og{SFRi) = log(5Fi?„(uv))+^(Auv) X 
, r 0.797 log(SFRi) -1-3.834 
2:88 



(7) 
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where X(365nm) = 2.061, X(280nm) = 2.512 
and X(150nm) = 3.574. Recent results from 
Adelberger fc Steidel (2000) suggest, however, 
that observed (attenuated) UV luminosity at short 
wavelengths (A « 160 nm) is a very poor indica- 
tor of SFR. They find that while galaxies with 
high SFRs (and intrinsic UV luminosities) show 
increased levels of obscuration (see Figure 11 of 
Adelberger & Steidel 2000) consistent with our 
primary assumption, the obscuration is such as 
to attenuate the observed UV luminosities to a 
remarkably constant level for a broad range of ex- 
tinctions (sec Figure 17 of Adelberger & Steidel 
2000). This results in a situation where Equation 7 
effectively becomes degenerate, and a knowledge 
of only the observed UV luminosity is insuffi- 
cient to regain information about the extent of 
the obscuration and the corresponding intrinsic 
luminosity or SFR. In contrast with this result. 
Figure 3 from Sullivan et al. (2001), comparing 
SFRuv derived from observations at A w 200 nm 
to SFRi.4GHz, docs not show the expected de- 
generacy, suggesting that perhaps the effect is 
stronger at shorter wavelengths. As we later show, 
we also find no such degeneracy in U-band de- 
rived SFRs, supporting this suggestion. The above 
method is thus explored below with regard to UV 
luminosities using the relations already derived, 
but placing emphasis on wavelengths longer than 
160 nm which appear less likely to show this de- 
generate effect. 

4. Corrections to local SFR estimates 

The relationship between (obscured) Ha- and 
UV-based measures of SFR and measurements 
of intrinsic SFR (for which we use SFRfir or 
SFRi.4GHz) are given by Equations 5 and 7. These 
relations are compared with observational mea- 
surements for a sample of local galaxies in Fig- 
ures 2 and 3, using data selected from the compila- 
tion of Cram et al. (1998). Figure 3 emphasises the 
validity of the radio-FIR correlation for this sam- 
ple, and justifies the strong emphasis on 1.4 GHz 
based estimates of SFR (and subsequently p*) in 
our analysis. The "knee" of the curves shown in 
Figures 2 and 3 marks the SFR at which the de- 
rived colour excess from Equation 3 goes to zero 
from higher values at higher SFRs. For SFRs lower 
than this value we assume there is no attenuation 
of the observed luminosity (and hence SFR). 



The SFR-dependent reddening formulation 
seems to account well for the observed discrep- 
ancy between SFRna and SFRfir or SFRi.4ghz 
(Figures 2a and 3a). It is also clear that, while 
this prescription may be useful for examining the 
trends over large samples, the extent of the intrin- 
sic scatter in the relevant correlations seen in the 
observations, is quite large, up to two orders of 
magnitude in some cases. Any detailed analysis 
of individual galaxies should thus be treated with 
caution, and quoted uncertainties in any derived 
results should reflect this. 

The attenuation of the UV-derived SFRs, how- 
ever, appears to be less effectively modelled. As 
shown by the dot-dashed lines in Figures 2(b) and 
3(b) the attenuation implied at U-band by Equa- 
tion 7 is insufficient to reproduce the observed 
SFRuv- The observations of local galaxies imply 
even greater levels of obscuration, and the degen- 
eracy in observed UV luminosities with respect to 
intrinsic UV luminosity found by Adelberger & 
Steidel (2000) possibly suggest greater levels still. 
The trend seen by Adelberger & Steidel (2000) 
could be illustrated in Figures 2(b) and 3(b) by 
a horizontal line at SFRuv ~ 1 Afoyr"^! for in- 
trinsic SFRs > lM0yr~^. This would intercept 
the solid line in these Figures at a value of in- 
trinsic SFR « 5M0yr~^, (less for the longer UV 
wavelength models) and the implied obscuration 
would be greater than given by Equation 7 (for 
A = 150 nm) for objects with larger intrinsic SFRs. 
Again, knowledge of only the observed UV lumi- 
nosity would be insufficient to establish the extent 
of the obscuration in the degenerate case. 

To establish whether or not the trend seen in 
the local U-band measurements is consistent with 
the results of Adelberger & Steidel (2000), we have 
applied several robust regression algorithms (Isobe 
et al. 1990) to the data of Figures 2(b) and 3(b). 
We find positive, non-zero slopes in all cases, even 
when considering only data having values of in- 
trinsic SFR > 1 Mq yr-i (the "flattest" portion of 
these diagrams). Our best estimate of the slope, 
from the ordinary least squares bisector (Isobe et 
al. 1990), is 0.7 ± 0.07 (for logarithmic values of 
SFR), with a Pearson correlation coefficient of 0.5. 
This may be compared with an approximate slope 
of 0.68 (for 0.1 <SFRfir < 10) for the relation 
predicted by Equation 7 at 150 nm {not U-band). 
This result will be made use of in Section 6. While 



5 



it is possible that some incompleteness in the data 
has artificially skewed the slope to positive val- 
ues, it is also possible that there is a real trend 
in the U-band data which is absent in the shorter 
wavelength (160 nm) data. This suggestion is not 
unreasonable, as the level of obscuration expected 
at U-band (365 nni) is less than at 160 nm for a 
given intrinsic luminosity, or SFR. It is also con- 
sistent with the non-degenerate trend seen in Fig- 
ure 3 of Sullivan ct al. (2001) based on 200 nm 
observations. If this is the case, then it is possi- 
ble that a method similar to the one we analyse 
here may be useful for correcting UV observations 
at wavelengths between U-band and 160 nm. In- 
deed, it should be emphasised that the empirical 
trends presented by the UV data (both in Fig- 
ures 2(b) and 3(b), as well as Figure 3 of Sullivan 
et al. 2001, and in Adelberger & Steidel 2000) are 
still evidence of a form of SFR-dcpendent obscu- 
ration, albeit one implying even greater levels of 
attenuation than predicted by Equation 7. 

5. Corrections to global SFR density 

We have so far considered the effects of a 
SFR-dependent obscuration on samples of galaxies 
taken from the local universe. Now, building on 
the conclusions of the previous Section, we con- 
sider how this might affect the evolution of the 
global SFR density, given the simplified assump- 
tion of no evolution in the properties or effects of 
dust. In light of the shortfall in the level of atten- 
uation predicted by Equation 7 we will treat the 
magnitude of the predicted corrections to as 
lower limits to the required correction, and exam- 
ine the resulting implications. 

In applying the SFR-dependent reddening cor- 
rection to samples at high redshifts, we need to 
assume that the locally-derived trend between ob- 
scuration and SFR is still valid. Obviously this 
may not be the case, particularly since our un- 
derstanding of galaxy evolution implies that the 
metallicity of the inter-stellar medium increases 
with time. For example, Steidel et al. (1999) find 
a dearth of high-redshift galaxies with E{B — V) > 
0.3. Extremely red galaxies, a class of objects de- 
fined hy R — K > 6, may actually be a popula- 
tion of dusty star-forming galaxies at high red- 
shifts (e.g. Cimatti ct al. 1999; Andrcani ct al. 
2000), however, so the extent of obscuration at 



high redshifts is still unclear (but sec also the dis- 
cussion of such objects in Adelberger & Steidel 
2000). In the absence of more extensive evolution- 
ary constraints, we begin with this simple assump- 
tion and anticipate that (1) evolutionary effects 
may be able to be incorporated as more data be- 
comes available, and (2) remaining discrepancies 
will serve to direct attention to those wavelengths 
and redshifts deserving of further investigation. 

5.1. Applying the correction 

We present a compilation of uncorrected global 
SFR density (p*) measurements in Figure 4. This 
emphasises the magnitude of the attenuation cor- 
rections required to reconcile the various estimates 
at different wavelengths. In addition to the UV 
surveys (Connolly et al. 1997: Trcycr ct al. 1998; 
Steidel et al. 1999; Sullivan et al. 2000) and Ha 
surveys (Gallego et al. 1995; Tresse & Maddox 
1998; Yan et al. 1999; Hopkins et al. 2000), we 
also show results from FIR observations (Rowan- 
Robinson et al. 1997; Flores et al. 1999) and 
1.4 GHz radio surveys (Haarsma ct al. 2000; Ser- 
jeant et al. 2000; Condon 1989). All measurements 
have been converted to our assumed cosmology 
and SFR calibrations for consistency. 

To correct existing UV and Ha estimates of p* 
using the SFR-dependent reddening formulation 
described above requires a knowledge not just of p* 
itself, but of the appropriate Ha or UV luminosity 
function, since the corrections to be applied are 
SFR (and hence luminosity) dependent. We have 
initially considered our derived relations, applying 
the corrections derived from Equations 5 and 7 
for a number of published Ha and UV surveys, 
sampling a broad redshift range. The results are 
shown in Figure 5. 

The local 1.4 GHz measurement from Serjeant 
et al. (2000) of p|.4 = O.O31±O.OO7M0yr-iMpc-3 
(after converting to Hq = 75 and using our as- 
sumed SFRi.4GHz calibration), derived from a 
measurement of the local 1.4 GHz LF of Revised 
Shapley-Ames spiral galaxies, is consistent with 
the value of 4 = 0.036 M0yr"^Mpc"^ obtained 
by integrating the local 1.4 GHz LF for spiral 
galaxies from Condon (1989). 

Also shown in Figure 5 are two models for the 
evolution of p* . The solid line is derived from a fit 
to the evolving 1.4 GHz LF for star- forming galax- 
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ies, measured by Haarsma et al. (2000) and based 
on the local LF of Condon (1989), and is anchored 
at 2; = by that point. The dashed line is a model 
for the FIR luminosity density derived by Gispert 
et al. (2000) from measurements of the cosmic far- 
infrared background. This model has been con- 
verted to SFR density using the SFR calibration 
of FIR luminosity given by Kennicutt (1998). 

5.2. Discussion 

The effectiveness of correcting estimates of p* 
by assuming an SFR-dependent obscuration is dis- 
cussed separately below for three broad redshift 
regimes. 

5.2.1. < ^ < 0.3 

For redshifts out to 2 « 0.3 applying the SFR- 
dependent correction to (1) the uncorrected UV 
LFs of Sullivan et al. (2000) and Treyer et al. 
(1998), and (2) the N[ii] and aperture corrected 
(but not extinction corrected) Ha LF of Tresse 
& Maddox (1998) all give consistent results. The 
magnitude of the correction for p^y in this red- 
shift range is a factor of about 3. While this cor- 
rection may need to be considered a lower limit the 
true value is probably not much greater than this, 
since both corrected pyy estimates and the cor- 
rected pjj^ estimate are consistent with the trend 
in the radio-based p* measurements spanning this 
redshift range (sec also Section 6). 

The one remaining discrepancy is a factor of 
2 — 3 difference between the 1.4 GHz and Ha based 
estimates of p* ai z k, 0. The results of Gallego et 
al. (1995) are shown in Figure 5 with the redden- 
ing correction originally applied by those authors. 
Mobashcr et al. (1999) have also derived a local Ha 
SFR density which is consistent with this measure- 
ment. The SFR-dependent formulation results in 
a correction similar to the k, 1 magnitude often 
assumed for this local population, and is unlikely 
to increase by more than a few tens of percent the 
value of Gallego et al. (1995). (At higher redshifts 
the evolving Ha LF provides greater numbers of 
higher luminosity systems, which results in larger 
corrections to /Ohq-) 

Serjeant ct al. (2000) address this discrepancy 
by suggesting that a significant fraction of star 
formation occurs in the cores of giant molecular 
clouds, which would imply much greater obscura- 



tion than derived from a Balmcr decrement red- 
dening correction assuming a simple dust screen. 
Regardless of the mechanism, however, the obscu- 
ration must be such that it explains the observed 
trends of Figures 2 and 3. Hence, at least in 
the local universe, the correction cannot be much 
greater than that which we examine in this anal- 
ysis. At higher redshifts, perhaps, the amount of 
attenuation could be greater, although the trend 
of Figure 2(a) is sufficient to explain the discrep- 
ancy in at 2; « 1 (see Section 5.2.2). 

An alternative explanation for this discrepancy 
questions the calibration of SFRi.4ghz- This is ul- 
timately based on the calibration of non-thermal 
radio luminosity, Lat, to supernova rate, vsn, in 
our own galaxy (Condon 1992). If the escape 
of cosmic ray electrons from our galaxy is sig- 
nificant, lowering the observed ijv for a given 
SFR, then the derived calibration constant will 
be larger than had we seen all the radio emis- 
sion. Hence, in external systems where we may see 
more, even most, of the radio emission produced 
by the star-formation, the derived SFR will be too 
great. Condon (1992) argues against this scenario, 
however, stating that significant variations in the 
ratio Ln/vsn would violate the observed radio- 
FIR correlation. 

A third possibility is suggested by preliminary 
estimates of the local Ha luminosity density from 
the KPNO International Spectroscopic Survey 
(KISS, Gronwall et al. 1997, 1999; Salzer et al. 
2000). This survey finds a local Ha luminosity 
density somewhat higher than that measured by 
Gallego et al. (1995). 

5.2.2. 0.3 < 2 < 2.0 

Within this redshift regime there is a discrep- 
ancy between estimates of Ppm from different 
analyses. For 0.3 < z < 1.0, the FIR-based values 
of Rowan- Robinson et al. (1997) are comparable to 
the radio-based results of Haarsma ct al. (2000), 
but those of Flores et al. (1999) lie a factor of 
« 1.5 — 2 lower (Figure 4). Possible sources for this 
discrepancy have been discussed by Flores et al. 
(1999), who emphasise that their analysis covers a 
much larger area than that of Rowan- Robinson et 
al. (1997), and uses different methods for fitting 
galaxy spectral energy distributions. 

From 0.6 < 2; < 2.0, the application of the 
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SFR-dcpcndcnt reddening to Ha LFs (Yan et al. 
1999; Hopkins et al. 2000) results in SFR densi- 
ties which are consistent with the observational 
radio-derived values, and the FIR-derived values 
of Rowan- Robinson et al. (1997). This is highly 
encouraging, and suggests that the form of SFR- 
dependent reddening described by Equation 5 may 
be applicable as far out as z « 2. It also implies 
little evolution in the form and extent of obscu- 
ration, at least in a global, population-averaged 
sense. Obviously details for individual galaxies 
will vary significantly, as evidenced by the large 
scatter in Figures 2 and 3. The integration of the 
luminosity function to derive estimates of p* has 
the effect of averaging over this scatter, minimising 
its effect. The resulting integral measure is thus 
more robust to the variations in individual galax- 
ies, and reflects instead the general trend such as 
that modelled by Equation 5. A further observa- 
tional estimate of at -^^ ~ 0.9 (Glazebrook et al. 
1999, not shown in Figures 4-6) may also be con- 
sistent with these results, as they comment that 
their value is consistent with the "high" estimate 
of Rowan-Robinson et al. (1997) if they apply the 
Calzetti prescription for dust attenuation. 

The same treatment with the UV LFs of Con- 
nolly et al. (1997) give a factor of w 4 corrections 
to Puv The corrected still lies a factor of 
« 1.5 — 2 lower than /OpjR, Pi.4ghz cor- 
rected (Figure 5). This is not unexpected, 
and is consistent with the interpretation of Equa- 
tion 7 as providing lower limits for the corrected 

It is possible, on the other hand, that lower 
mctallicities and obscurations at higher rcdshifts 
for galaxies of a given luminosity imply some evo- 
lution in the calibration of SFRyv- Such an effect 
would be in the sense that a given observed UV 
luminosity at higher redshifts implies lower lev- 
els of SFR. This would give a calibration coeffi- 
cient which decreases with redshift, reducing the 
observed Puyi and the correction for any obscu- 
ration would also decrease with redshift. If this 
effect is large enough, it could exacerbate the dis- 
crepancy between pyy and the 1.4 GHz and FIR 
estimates. This extreme scenario would be incon- 
sistent with the results of the Ha corrections in 
this redshift regime, which suggest that any such 
evolution in the level of obscuration is probably 
small. 



5.2.3. z > 2.0 

At the high redshift end, z « 3 - 4, the UV LFs 
of Steidel et al. (1999) have also been re-examined 
assuming an SFR-dependent reddening. Despite 
the large uncertainties associated with extrapolat- 
ing Equation 7 to high redshifts, and the expec- 
tation that in its present form the correction is 
likely to be underestimated, it is encouraging to 
see that the results lie midway between the pre- 
dictions of the two models shown. The corrected 
estimates using Equation 7 are bracketed by the 
predictions from the model evolving radio LF and 
the FIR background model. They are also con- 
sistent with the lower limit derived from sub-mm 
SCUBA observations of the Hubble Deep Field 
(Hughes et al. 1998), p* > O.ieMgyr-^Mpc-^ in 
our cosmology. The correction derived from Equa- 
tion 7 approaches a factor of 10 at z ~ 4, and 
may still need to be considered a lower limit to 
the true correction given the results of Adelberger 
& Steidel (2000). The extremely large uncer- 
tainties associated with extrapolating the locally- 
derived SFR-dependent reddening models to such 
high redshifts, though, dictate the use of caution 
when interpreting these results. It is still inter- 
esting to note that if these corrections are indeed 
underestimates and the true values need to be 
higher, they would lie closer to the extrapolations 
from the 1.4 GHz luminosity fimctions estimated 
by Haarsma et al. (2000), almost an order of mag- 
nitude above the estimates from the FIR back- 
ground. 

6. An empirical investigation 

The inadequacy of Equation 7 in fully explain- 
ing the observed levels of obscuration prompts two 
questions. The obvious one relates to the origin 
of the obscuration, and is discussed briefly below. 
The second is: "Independent of the mechanism 
producing the obscuration, is the local empirical 
relation between SFRuv and intrinsic SFR sufB- 
cient to explain the discrepancies in p* at higher 
redshifts?" We examine this by applying the em- 
pirical relationship for local galaxies observed in 
Figures 2(b) and 3(b) to measurements at all 
redshifts. 

For this simple analysis we treat the relation- 
ship between SFRuv and SFRfir (or SFRi.4ghz) 
described by the 0.15 /xm attenuation (the solid 
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curve in these Figures) as a mathematically con- 
venient description of this empirical relation. We 
also assume this relationship is valid at shorter 
wavelengths (where the attenuation may be still 
greater, given the results of Adelberger & Steidel 
2000), but particularly at 0.28 /zm, where the mea- 
surements from Connolly et al. (1997) were made. 
Applying the SFR-dependent correction to all the 
Puv points based on this empirical result gives the 
diagram shown in Figure 6. Good agreement is 
now seen between all UV, Ha and radio-derived 
values for p* . The effectiveness of this application 
is dominated by the effective correction to the data 
of Connolly et al. (1997), since the corrections at 
low-redshift are not modified greatly (due to the 
LF being dominated by lower- luminosity sources) . 
The primary conclusion here is that the empirical 
relation between local U-band estimates of SFRuv 
and intrinsic SFR is sufficient to account for ob- 
scuration in 280 nm observations at z f» 1.5. 

Two possible scenarios to explain the short- 
comings of Equation 7 are discussed in this light. 
First, dust in star forming galaxies may be char- 
acterised by a reddening curve even greyer than 
that of Calzetti et al. (2000). The situation may 
be explained if the attenuation at wavelengths be- 
tween 0.15 < A < 0.4 /im is comparable to that 
observed locally in U-band. One model produc- 
ing this effect may be a heavily dust-enshrouded 
starburst, optically thick below 0.4 /um, but having 
"holes" which allow some portion of UV light to es- 
cape imattenuated. This would result in a uniform 
"blocking" of all emission below the optically thick 
cutoff, and produce the required constant attenu- 
ation as a function of (UV) wavelength. The ef- 
fective reddening curve for this scenario may, how- 
ever, not be consistent with the results of the FIR 
dust luminosity analysis presented in Calzetti et 
al. (2000). Alternatively, the dust attenuation re- 
mains as given by Calzetti et al. (2000), but there 
is some additional process which needs to be in- 
voked to explain the shortfall in the U-band de- 
rived SFRs seen in Figures 2(b) and 3(b). This 
may be the case, for example, if the relative contri- 
bution of the old stellar population to the U-band 
luminosity is greater at low SFR. (Luminosities 
at shorter UV wavelengths would have less "con- 
tamination" from the old stars, so would be af- 
fected less or not at all in this scenario.) This 
would imply that the U-band derived SFRs in low- 



SFR systems are ower-estimated, and need to be 
revised downward. Then, combined with a suit- 
able revision of the U-band/SFR calibration, this 
would account for the discrepancies. There are 
still problems with this suggestion as well, since 
low-SFR systems also tend to be low-luminosity 
systems (e.g. Cram et al. 1998, their Figure 3), 
and this implies the relative contribution from the 
old stellar population may not change much with 
SFR. 

In addition, whatever mechanism is invoked 
needs to be able to reproduce the degenerate rela- 
tionship between observed UV luminosity and ob- 
scuration at shorter wavelengths found by Adel- 
berger & Steidel (2000). More investigation of 
the result described here is obviously still re- 
quired. Although the application of the locally 
observed empirical SFRuv to SFRfir relation to 
shorter wavelengths and higher redshifts may ac- 
count for the discrepancies it does not explain 
their physical origin. It does emphasise that the 
trend observed locally in U-band is consistent with 
the discrepancies at much higher redshifts and at 
shorter wavelengths. This is consistent with the 
result for the Ha correction suggesting little net 
evolution, if any, in the level of obscuration at least 
out to ^; « 2. 

7. Conclusions 

We have modelled an SFR-dependent attenu- 
ation by dust, characterised by the Calzetti red- 
dening curve. Corrections derived from this model 
have been applied to a large sample of local galax- 
ies, spanning a wide range in intrinsic SFR, to 
examine the effects on local Ha- and UV-derived 
SFR estimates. Subsequently, measurements of 
p* for a broad range of redshift have also been 
corrected using this model. It is clear from the 
observational data investigated here, as well as 
the smaller samples of Wang & Hcckman (1996) 
and Calzetti et al. (1995), that some form of 
SFR-dependent extinction is implied. The inves- 
tigated prescription appears sufficient to explain 
the general trend between SFRhq and SFRi.4ghz 
or SFRfir. It also accounts for the discrepancies 
between estimates of p'^^ and p\ 4ghz or PpiR out 
to » 2, with the exclusion of the local p'^^ es- 
timate of Gallego et al. (1995). This is strongly 
suggestive of little or no evolution in the extent 



9 



or form of the obscuration at rcst-framc Ha wave- 
lengths. As a result, this may serve as a useful tool 
in making corrections for obscuration in estimates 
of SFRho or in the absence of more direct 
methods (such as for the NICMOS grism surveys 
of Yan et al. (1999) and Hopkins et al. (2000), 
and the recent SOFT/ISAAC results of Moorwood 
et al. (2000), which lack Balmer decrement infor- 
mation). 

The examined prescription for correcting UV- 
based estimates of SFR and p* was less effective. 
Although the predicted corrections reduce the 
level of observed discrepancy with obscuration- 
free estimators, they do not fully account for the 
observed levels of attenuation in local U-band es- 
timates of SFRuv- There may also be concerns 
that a degenerate relation (at shorter wavelengths) 
would prevent application of the iterative method 
used. The use of this method to estimate lower 
limits for the corrections to pjjy is more encour- 
aging, though, particularly for z < 0.3 where the 
magnitude of the required correction is unlikely 
to increase by much. The situation becomes more 
complex at higher redshifts, and here evolutionary 
effects, as well as a more complete understanding 
of the obscuration in galaxies at UV- wavelengths, 
may need to be incorporated. 

The empirical relationship observed at U-band 
for local galaxies between SFRuv and SFRfir or 
SFRi.4GHz was also examined. It was used to cor- 
rect high-redshift estimates of Puv- resulting in 
values consistent with p*.4ghz ^^'^ Pfir ^'^^ UV 
wavelengths longer than 160 nm. This may sup- 
port the result from the Ha analysis, suggesting 
minimal evolution in the extent of dust obscura- 
tion. The mechanisms producing the obscuration 
are still in question, and prompt further investiga- 
tion of dust models and the characteristics of dust 
on UV-wavelength radiation. 
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Fig. 1. The colour excess E{B — V) as a function of FIR luminosity. Solid triangles: data from Wang & 
Heckman (1996, their Figure 10 with Balmer decrement converted to colour excess); Open circles: starburst 
and blue compact galaxies from Calzetti et al. (1995). The solid line shows the relation between colour 
excess and FIR luminosity derived from combining Equations 1 and 2. 
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SFRp,jj (Mq yr ') SFR^^^ (M^ yr ') 

Fig. 2. — (a) SFRhq, and (b) SFRuv, both uncorrected for extinction, compared to SFRfir. The SFRfir 
values here come from measurements of the 60 /xm luminosity, and the SFRuv values come from U-band 
(365 nm) measurements. Both of these are converted to SFR using the calibrations of Cram et al. (1998), 
after adjusting by a factor of 5.5 to account for stellar masses < 5Mq. Both calibrations are consistent 
with those of Kcnnicutt (1998). The solid line shows the SFR-dependent attenuation, calculated using the 
prescription described in the text, for (a) Ha and (b) UV(A = 0.15 y^m). The dashed line shows a one-to- 
one relationship. The dot-dashed and dotted lines in (b) are for the attenuations valid at A = 0.365 fim 
(U-band) and 0.28 /im respectively. It can easily be seen that the U-band attenuation in this formulation 
is not siifficient to reproduce the observed trend. The different symbols mark different sources for the data 
(see also Cram et al. 1998): KK83: Kcnnicutt & Kent (1983); KKHHR87: Kennicutt et al. (1987); K92: 
Kennicutt (1992); R90: Romanishin (1990). 
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Fig. 4. — SFR density (p*) as a function of redshift. This diagram is a compilation of SFR densities derived 
from various existing sources, with no reddening corrections applied to any data. The reddening-corrected 
measurement of Gallcgo ct al. (1995) has been artificially "reddened" by 1 magnitude for comparison with 
other uncorrected data in this diagram (solid circle). References in diagram are as follows, along with the 
origin of the p* estimate: HaaOO: Haarsma ct al. (2000) (1.4 GHz); HopOO: Hopkins et al. (2000) (Ha); 
SerOO: Serjeant et al. (2000) (1.4 GHz); SulOO: Sullivan et al. (2000) (UV); Flo99: Florcs ct al. (1999) (FIR); 
Ste99: Steidel et al. (1999) (UV); Yan99: Yan et al. (1999) (Ha); Tre98: Trcyer ct al. (1998) (UV); TM98: 
Tressc & Maddox (1998) (Hq); Con97: Connolly et al. (1997) (UV); RR97: Rowan- Robinson et al. (1997) 
(FIR); Gal95: Gallego et al. (1995) (Ha); Con89: Condon (1989) (1.4 GHz). 
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Fig. 5. — SFR density (p*) as a function of redshift. This diagram shows the same data as in Figure 4 
with the Ha and UV based measurements corrected using the SFR-dependent reddening prescription. The 
FIR measurements of Florcs ct al. (1999) have been omitted in this diagram for clarity. The solid curve 
comes from the evolving 1.4 GHz radio luminosity function for star- forming galaxies derived by Haarsma et 
al. (2000). The dashed curve is a model derived from the cosmic far- infrared background by Gispert et al. 
(2000). References are as for Figure 4. 
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Fig. 6. — SFR density (p*) as a function of redshift. This diagram is identical to Figure 5 except that the 
values have all been corrected using the SFR-dependent reddening with A;(A) = fc(0.15 /xm), to examine the 
effects of applying the locally observed empirical relationship between SFRuv and SFRfir to all redshifts. 
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